Effective temperature of active matter 
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We follow the dynamics of an ensemble of interacting self-propelled particles in contact with an 
equilibrated thermal bath. We find that the fluctuation-dissipation relation allows for the definition 
of an effective temperature that is compatible with the results obtained using a tracer particle as a 
thermometer. The effective temperature takes a higher value than the temperature of the bath when 
the effect of the motors is not correlated with the structural rearrangements they induce. Instead, it 
takes a slightly lower value than the environmental temperature when susceptible motors are used. 



Active matter is driven out of equilibrium by internal 
or external energy sources. Its constituents absorb energy 
from their environment or from internal fuel tanks and 
dissipate it by carrying out internal movements that lead 
to translational or rotational motion. A typical example 
is given by the eukaryotic cell cytoskeleton network [l| 
that is made of long polar filaments with attached mo- 
tors. The latter exert stresses which deform the former 
and regulate the network dynamics. These structures 
exhibit a rich variety of viscoelastic properties; by rear- 
ranging their structure they change from plastic/fluid to 
elastic/solid phases and vice versa. Another celebrated 
example of active matter are self-propelled particle as- 
semblies in bacterial colonies [2J. 

In this Letter we address the following fundamental 
question: which thermodynamic concepts, if any, can be 
applied to active matter? We focus on the definition of 
an effective temperature, T e ff } through the comparison of 
induced and spontaneous dynamic fluctuations [3]. The 
notion of such a T e f f was introduced in the study of (pas- 
sive) glassy systems, that is to say, macroscopic objects 
with a sluggish relaxation that takes them to an asymp- 
totic small-entropy production regime which is, though, 
still far from equilibrium. In these systems T e fj is a good 
thermodynamic concept. It admits a microcanonical-like 
definition in which only asymptotically dynamically ac- 
cessible states contribute to the entropy [4[ • T e f f can be 
measured by a fine-tuned thermometer, it satisfies a zero- 
th law in the sense that it takes the same value for all ob- 
servables evolving in the same time-scale, and heat flows 
in the direction of negative effective temperature differ- 
ences Q. Moreover, it appears in the modification of 
fluctuation theorems when applied to systems that can- 
not equilibrate when let freely evolve [5[ . These features 
were confirmed with numerical simulations of a variety 
of realistic glassy models Q and a number of groups are 
testing these ideas experimentally on colloidal suspen- 
sions [7| . Theory suggests that the effective temperature 
should also be well-behaved in weakly driven systems - 
in a small entropy production regime 3]. This fact was 
put to the test in gently sheared super-cooled liquids and 
glasses Q , vibrated granular matter Q , and moving vor- 



tex phases in superconductors Q . In biologically inspired 
problems the relevance of T e j f was already stressed in 
studies of gene networ ks J 1 01] and to reveal the active 
process in hair bundles [111 ]. 

The role played by Tejf in the stability of dynamic 
phases of motorized particle systems was stressed by Shen 
and Wolynes [l2| with a variational analysis of the mas- 
ter equation of a similar model to the one we study here. 
In order to investigate the existence and properties of 
T e ff in active matter we study a schematic statistical 
physical model: a system of i = 1, . . . , N massive, in- 
teracting and motorized spherical particles confined to a 
fixed volume, V, in 3d. Spherical particles are indeed too 
simple for real applications in the context of cytoskele- 
ton studies but we choose them to highlight the effect of 
the non-equilibrium drive. Similar behaviour is expected 
for other geometries of the constituents. Particles are 
in contact with a thermal environment that is described 
by a random noise and a viscous drag. The particles' 
velocities Vj thus evolve with the Langevin equation: 



miVi = -£miVi + ff + f, 



(1) 



the frictional force and £ the fric- 
r\i is a Gaussian white noise represent- 



Here, — ^rriiVi is 
tion coefficient, 
ing thermal agitation with zero mean, and correlations 
(Vi(t)r)i(t')) = 2£miTd(t-t') (we set k B = 1 henceforth). 
The systematic mechanical conservative force on particle 



i due to all others is ff = Y?j=i tfj = _ ^i^( r v) 
with the 2n-n Lennard-Jones potential [13] U(r 



4e[(a/rij) 2n — (cr/r^)"]. is the interparticle distance. 
The short-range repulsion, parametrized by a, prevents 
the overlapping of the particles and it is then a (soft) 
measure of the particle diameter. The energy parame- 
ter, e, describes the depth of the attraction. The mid- 
range attraction, rMin = cr2~ 1 / n , mimics the effective 
cross-linking by linker proteins between elements in the 
cytoskeleton. By choosing n — 18 we get an adhesive 
soft-sphere potential with rMin <C o as often used to de- 
scribe biomolecular assemblies 12 , [13] • Finally, motors 
apply a force f 4 M on each particle following a stochastic 
process that mimics the realistic chemical process. For 
simplicity, we use here a time series of isotropic kicks. 
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FIG. 2: (Colour online) Main panel: the incoherent inter- 
FIG. 1: (Colour online) The static structure factor, S(Q), mediate scattering function at Q = Qntax, for motor force 
for a system in equilibrium and driven out of equilibrium by intensities given in the key. Inset: the dependence of the a 
motors with / = 50%F. relaxation time on the strength of the motor forces with an 

exponential fit r a ~ Aexp -(/ M //°) with f° ~ 0.90. 



During r steps of the molecular dynamics trajectory in- 
dependent forces are applied to a fraction of randomly 
chosen particles. The strength of the force exerted on 
each particle is chosen to be a fraction of the mean me- 
chanical force acting on the equivalent passive system, 
F = jj YljLi T w0 types of motors are considered: 
adamant and susceptible [121 ] . The former act in such a 
way that the direction of the force is chosen at random 
isotropically: their action is independent of the struc- 
tural rearrangements they induce. The subset of pro- 
pelled particles and the directions of the applied forces 
change at each power stroke of duration r. Susceptible 
motors slow down when they drive the system uphill in 
the (free)-energy landscape. In both cases motors are 
permanently applied. We neglect hydrodynamic effects. 

We integrated Eqs. |T]) numerically using Ermak's al- 
gorithm [l5|, LL6( . We henceforth use standard Lennard- 
Jones units: length, energy and time-units are a, e and 
(mcr 2 /e)2, respectively [161 ] . In all our simulations we 
used 100 independent configurations with 500 particles. 

The passive system (motors switched off) is well char- 
acterized and presents gas, liquid and crystalline phases 
depending on the bath temperature and particle den- 
sity [HI]. By tuning the intensity and direction of the 
non-conservative forces we drive the system near equilib- 
rium (weak perturbation) or far from equilibrium (strong 
perturbation) . We delay the description of the dynamic 
phase diagram to a more detailed publication. Here we 
focus on the driven dynamics of a system at conditions 
such that it is a liquid in equilibrium: T = 0.8 and p = 1. 

Let us first discuss the effect of adamant motors. 
We show data using a time-scale for the power strikes 
t = 5 x 10 3 , we apply the force to 10% of the parti- 
cles and the strength of the non-conservative force, / , 
is indicated in each figure. We start by comparing the 



structure of active and passive matter. In both cases 
the structure factor, S(Q) = (p(Q,t)p*(Q,t)), with 
p(Q, t) = N- 1 Y£=i e'^iM the Fourier transform of the 
instantaneous density, becomes stationary and isotropic 
after a short transient: S(Q,t) — > S(Q). In Fig. Q] we 
show that S(Q) for the passive liquid and the active sys- 
tem with f M — 50%F are practically identical, with a 
maximum located at Quax ~ 7.1 in both cases. 

The relaxational dynamics of super-cooled liquids and 
glasses is usually examined by monitoring the decay of 
the - possibly two-time dependent - incoherent (one- 
particle) intermediate scattering function F s (Q,t,t w ) — 
N' 1 J2^ =1 (e- l ^ ( - t+t ^- p '^ ) 0] with t w the waiting- 
time measured after preparation and t the delay time 
between the total and the waiting times. In the liquid 
and under the effect of the motors the systems reach an 
isotropic stationary regime: F s only depends on Q and 
t. Figure [2] displays F s as a function of time-delay in a 
log-linear scale for several force strengths ranging from 
zero (equilibrium limit) to 90%F. The decay is faster for 
increasing force strengths though remains relatively slow 
for all drives as demonstrated by the fact that F s decays 
to zero in a logarithmic time-scale. In the inset we show 
how the a-relaxation time, r Q = F s (Q,T a ) — e~ l , de- 
creases for increasing force strength in a linear-log scale, 
similarly to the shear-thinning phenomenon. The line 
is r a ~ Aexp ~(f M / f°). At each applied force we find 
r a oc Q~ 2 (not shown). 

The equilibrium fluctuation-dissipation theorem 
(FDT) states that the spontaneous fluctuations are re- 
lated to their associated induced fluctuations by a model- 
independent formula T\(t, t w ) = C(t — 0, t w ) — C(t, t w ). 
C is the connected correlation of two observables 0\ and 
02, measured at times t + t w and t w respectively, x is 
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FIG. 3: (Colour online) Parametric representation of the fluc- 
tuation dissipation relation for passive and active matter. The 
strength of the applied forces are given in the key. 



the linear response of the average of the observable 0\ 
measured at time t + t w to an infinitesimal perturbation 
that modifies the Hamiltonian as H — » H — hOi from 
t w to t: X (t,t w ) = (Of(t + t w ) - O x {t + t w ))/h. 
In interacting particle problems [8j it is cus- 



to 



use 



Oi(t w ) - jfYsiLi 



ternary 

2 {t w ) = 2Y / i = i^cos[Qr i (t w )], with the field e, = ±1 
with probability a half. For each system configuration 
we averaged over 100 field realizations. In Fig. [3] we 
show the fluctuation-dissipation relation for passive and 
active matter. The plots are parametric constructions 
of the integrated linear response x as a function of the 
correlation function using t as the parameter [3]. In the 
equilibrium case, FDT holds and minus the inverse slope 
of x(C) is the temperature of the thermal bath, T = 0.8 
in this case. In general, all curves join at (C — 1, x — 0), 
as imposed by normalization at t = 0. The 'initial' 
slope is determined by the equilibrium FDT [T^. At 
longer time-differences the curves progressively depart 
from the equilibrium result to reach other straight lines 
characterised by slopes that depend on the strength of 
the motor forces. A time (or correlation) dependent 
effective temperature is then defined as Q 



T eff (C) = -d X (C)/dC 



(2) 



Interestingly enough, soon after departure from the equi- 
librium form, the curves approach a new straight line 
from which one extracts a single value of the effective 
temperature - within numerical accuracy [l8j ]. 

The meaning of T e ff as a temperature depends on 
it verifying a number of conditions expected from such 
thermodynamic concept. In particular, a temperature 
should be measurable with a thermometer. A tracer par- 
ticle with a long internal time-scale (proportional to the 
square root of its mass) d, 0, Q acts as a thermometer 
that couples to the long time-delay structural rearrange- 




FIG. 4: (Colour online) Main panel: probability distribution 
function of the velocity v of different tracers with masses mt r 
given in the key. The motor force intensity is / = 60%f\ 
The full lines are fits to the Gaussian form (|3}. Inset: T e ff 
as function of mt r as obtained from the Gaussian fits to the 
data-points in the main panel. 



ments (and not to the fast vibrations, t ~ 0, that yield 
the ambient temperature, see Fig. [3] and 17\). We then 
coupled a tracer particle with mass mt r to the active 
matter via the same Lennard- Jones potential, so that we 
do not modify the structure factor of the fluid. (The 
tracer does not couple to the thermal bath.) In Fig. 2] we 
show the tracers' velocity distributions (for 10 indepen- 
dent tracers). The data points are superposed to 
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(3) 



with T e ff the only fitting parameter. The reason why a 
Maxwellian distribution applies is that the tracer behaves 
as a normal system immersed in an environment made 
of active matter. The inset in Fig. [4] displays the depen- 
dence of T e f f on mtr for one value of the energy pumping 
force. A very light tracer basically follows the very fast 

high frequency - dynamics of its environment and thus 
measures the bath temperature (lower horizontal dashed 
line). A heavier tracer feels the slower - low frequency 
- structural relaxation and measures a higher tempera- 
ture. A sufficiently heavy tracer only follows the slow 
dynamics of the sample and therefore measures the ac- 
tual effective temperature (upper horizontal dashed line) , 
T eff (f M = 60%F) - 1.23. 

Finally, in Fig. Owe compare the values of T e ff de- 
rived from Eqs. ([2]) and ([3]). The two measurements 
yield consistent results, with T e ff increasing from T at 
zero pump to about 2T when the strength of the driving 
forces equals the one of the mechanical forces. We also 
found that T e f f decreases with decreasing r (not shown) . 

Susceptible motors 12] act differently: they slow down 
when driving the system uphill in the (free)-energy land- 
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FIG. 5: (Colour online) Effective temperature obtained from 
response-correlation and massive tracer calculations for dif- 
ferent motor force intensities. 



scape. It is natural to expect that the effective tempera- 
ture of such active matter would be lower than the one of 
the thermal bath. We tested the effect of a very simple- 
minded realization of these motors. For r time steps the 
motors attempt to act on a percentage of the particles. 
Their effect is accepted only if the move leads to a de- 
crease in potential energy. Otherwise, the effect of the 
motors is neglected and free-relaxation takes place, fol- 
lowed by a new attemp to drive the system. It is clear 
that the number of accepted moves under the effect of 
such motors decreases with time and, in practice, the 
system is not sufficiently driven out of equilibrium. Fig- 
ure [6] shows a tiny effect, T e // S T. A more detailed 
study of susceptible motors will be presented elsewhere. 

Summarizing, we showed that the notion of an effective 
temperature, defined by studying the deviations from the 
equilibrium FDT out of equilibrium can be applied to ac- 
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FIG. 6: (Colour online) The parametric plot of linear response 
against correlation for active matter under the effect of sus- 
ceptible motors. 



tive matter. We demonstrated that the value of the effec- 
tive temperature measured with a thermometer realized 
by a tracer particle coincides - within numerical accuracy 
- with the one obtained from the fluctuation-dissipation 
relation for all values of the pumping forces. For adamant 
motors T e f / is larger or equal than the bath temperature 
and increases as a function of the pumping force. This 
finding is consistent with the intuitive idea that ascribes 
T e ff to the degree of additional agitation caused by the 
non-conservative forces. Susceptible motors have a very 
weak effect in driving the system out of equilibrium. For 
intermediate time-delays the fluctuation-dissipation rela- 
tion indicates T e f f < T but a much more detailed study 
is necessary to grasp the dependence on the interven- 
ing parameters such as the time-delay between kicks, the 
strength of the perturbing force, etc. in this case. Finally, 
it is a challenge to derive this results with hydrodynamic 
theories of the kind proposed in . 

LFC is a member of Institut Univcrsitaire de France. 
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